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Abstract 

In the framework of the quasipotential method we calculate the 
contributions of the kind a^lna, of three-photon exchange dia- 
grams to the energy spectrum of muonium v?S\ state in the leading 
order on rrie/m^. Analytical expression of obtained correction for arbi- 
trary principal quantum number is equal to (Za)^(me/?7i^)me(61n(2) — 
35/48) /n^. 
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The investigation of muonium and positronium fine structure represents 
one of the basic tests of quantum electrodynamics, which is sensitive to radia- 
tive corrections of higher order on a [|T|. Many papers 0, ^, Q are devoted 
to the calculation of different contributions to the fine structure of hydrogen- 
like atom energy levels. The interest to this problem remains unchanged 
[|, 1^, If. The progress, achieved in the last years during the calculation of 
logarithmic contributions of order In a in the positronium fine structure 
intervals {2^Si ^ l^Si, 2^Si ^ 2^Pj) |, [10|, [iTl, doesn't abohsh the necessity 



of calculation of higher order corrections 0{a^) |jT2|. The development of 
experimental methods, based on Doppler-free two-photon spectroscopy, al- 
lows the " large" structure intervals to be measured for the muonium and the 
positronium [O, |T1|, 113 1. 



. .3c. ^^ 1233607218, 9 ±10,7 Mi/z 
AEp^ (2 6i ^ 1 S,) - j ^233607216, 4 ± 3, 2 MHz 

AElZifSi ^ V^Si) = 2455527936 ± 120 ± 140 MHz. (2) 

The increase of the experimental accuracy of muonium fine structure 
interval measurements (just as for positronium), planned in the near future, 
makes very actual the calculation of radiative corrections of higher order on 
a. In this paper we have performed studies of the contributions of order 
O(a^lna), 0{a^) from three-photon exchange diagrams to the muonium 
fine structure. The contribution of these diagrams to the muonium hyperfine 
structure was obtained in [|16|. Our calculations are based on the Schrodinger- 



type local quasipotential equation ||T7 



^m{p) = / jT^^Vip, g, M)^M(g1, (3) 



\2iJ,R 2^r) J (27r) 

where IP' = El — m\ = E"^ — m|, iir = E1E2/M is relativistic reduced mass, 
M = El + E2 is the bound state mass, mi, m2 are the masses of the electron 
and the muon. As an initial approximation of quasipotential V{'p, M) for 
the bound state system (e~yU+) we choose the ordinary Coulomb potential. 



On the basis of equation (3) in [|T^ we have obtained some relativistic correc- 
tions ma^ in the positronium fine structure from the one-photon, two-photon 
interaction and the second order perturbation theory. There are six diagrams, 
shown on Fig.l, which determine the three-photon exchange interaction in 
the muonium. 

Let consider the first diagram of Fig.l. The corresponding amplitude 
already has the factor a^, which appears due to electromagnetic vertices and 
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Figure 1: Feynman diagrams of three-photon exchange interaction 
in the system (e~//+). 

Coulomb wave function. So, in the first stage of our calculations wc have 
neglected by the electron and muon vector momentum of relative motion in 
the initial and final states, taking into account that the necessary accuracy is 
already achieved. Then the first diagram amplitude (Fig.l) takes the form: 

^ jZaf r^^ r < ^^{gi - p' + mi)7i^(pi - p + mi)7f > 
47r^ J J {p^ — w"^ + ie) {p'^ — w'^ + ie) [{p — p'Y -\- ie\ 

< 72 (P2 + P + ^2)72 (g2 + + ^2)72 > 
D,ip)D,ip')D,{-p)D,{-p') 

where Df.,n{p) are the denominators of electron and muon propagators: 

D{±p)=p^ -w'^ ±2mp^ + ie, = -^)^ (5) 

and the corner brackets designate the averaging on Dirac bispinors; pi,p2 
are the four-momenta of particles in the initial state; qi^q2 are the particle 
four- momentum in the final state. As usually is, the factor Za emphasizes 
exchanging character of photon interaction between particles. The exchange 
photon propagators were taken in Feynman covariant gauge. As it famous. 
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using of Coulomb gauge is the most natural for exchange photons, because 
the Coulomb interaction dominates in the system (e"/!"*"). Nevertheless, the 
equivalence of Coulomb and Feynman gauges in the scattering approximation 
for the three-photon diagrams calculations was shown in |[T^. To construct 
the quasipotential of the system (e"/!"*") with L=0 and J=l, that corresponds 
to Ti'^ , let introduce the projector operator for initial and final states of the 



following kind |18, 19 



1 P + M . 

^ = 7= 77 ^5 (6) 

where P = pi + P2 = Qi + (I2 is, full four-momentum of the system, and 
is the vector of ^Si muonium polarization. Projecting the particles on ^5*1- 
state by means of (6), we avoid cumbersome matrix multiplication in the 
bispinor averages and immediately pass on to calculation of total trace in 
(4). As a result, the quasipotential of first diagram may be written in the 
form: 

3, ^ _{Zal r , r , F,{p,p') 

vrs J n ^ D,{p)D,{p')D,{p-p')D,ip)D,{p')D^{-p)D^{-p') 

(7) 

-D^(p) = p"^ -w"^ + ie, 

Fi {p, p) = fi2{p, p')ml + ^/iima, /12 = pp - ^rri\ - 2mipo - 2mipQ - 2poPo, 

(8) 

/ii(P)P') = 2mip''^ + pop'"^ + lOmipp' + 2popp' + 2pQpp'+ 

2 / 2 2 2 / 2 /2 / 

2mip + PqP + Qm-^^po + Gm^pQ + AttiiPq + Aniip^ — AitliPqPq. 

We kept in (7) only the terms proportional to m| and m2, taking in mind 
the determination of contribution to muonium fine structure in the leading 
order on parameter mi/ 1712- As will soon become evident, we can't restrict 
in F{p,p') only by terms ~ ml. The quasipotentials of the rest amplitudes 
of Fig.l may be constructed in a similar way. They differ from each other 
due to momentum dependence in muonic denominators and to the kind of 
functions fa (i=l,...,6). 

The parts of Fi{p,p'), proportional to m^, coincide in all six amplitudes. 
Let remark, that when substitute e — ^ 75 in projector operator (6) (^5*0 state), 
we obtain the same function fi2{p,p'), as for the ^Si muonium. This means, 
that the muonium hyperfine splitting appears as effect of higher order on 
mi/m2. Functions fn are equal to: 

/21 = —lOmip'^ — 5pop''^ + 10m ipp' + Ap^pp — Ap'^pp' + 2m + Ap'^p^ + 1 2pQm\ — 

(9) 
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-QtuIpq + Arriipl + SmipoPo ^ Smip'o^ + 4poPo, 
/31 = 2m\p'^ + 4pop'^ + IQrriipp' — ApQpp' + 'ip'^pp' — lOmip^ — Spgp^— (10) 

/41 = 2mip'^ +PoP'^ - 2mipp' + 4poPp' + 2poPp' - 10m ip^ - 8poP^ - 12mlpo+ 

(11) 

+6m^Po + 4miPo - 4miPoPo + ^miPg - Sp'qpI, 

/51 — lOmip'^ — 8pop'^ — 2mipp' + 2pQpp' + Ap'^pp' + 2mip^ +p'oP^ + Qmlpo— 

(12) 

/ei = — lOmip'^ — bpop''^ — 2mipp' — 4:pQpp' — -ip'^pp' — lOmip^ — hp'^p^— (13) 
—Qm\pQ — Qm\p'f^ — SmiPg — 4miPoPo ~ 8"^iPo • 

Integral function in (4) has simple poles on loop energies poiPo in electron, 
muon and photon propagators.. So, the most natural way of integration (4) 
consists in the calculation of integrals on po,Po ^^^^ step, using the 

method of residues. But such an approach of calculation leads, nevertheless, 
to rather complicated intermediate expressions, what makes highly question- 
able its subsequent analytical integration on spatial momenta p, ff . So, we 
have used different approach of integration in (4), connected with trans- 
formation of muonic denominators, accounting the necessary calculational 
accuracy on parameter m\/m2. Considering that the spatial momentum of 
muonic motion in the intermediate state \p\ < m2, we obtain: 

Di^ip) =p^ -nP + 2m2Po ~ 2m2 ipo — h ie j !^ 2m2(po + i^), (14) 

where the second approximate equality means that we neglect by the muon 
kinetic energy in the intermediate state. Doing so, we suppose that the in- 
tegration contour on variable po must be closed in the lower halfplane. Con- 
sidering the terms, proportional to in the numerators of all six diagrams 
(function fi2{p,p')), we have arrived to the need of expression transforma- 
tion, which includes the sum of muonic denominators (5). Using the second 
approximate equality from (14), we obtain: 

111 

f 



D^{-p)D^{-p') D^{-p)D^{p' - p) D^{-p')D^{p - p' 
1 11 



D^{p)D^{p - p') D^{p')D^{p' - p) D^{p')D^{p) 
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~ 2m2 2m2 ' 

In the energy spectrum the expression (15) will cause the corrections 
of order 0(q;^), which are canceled by the similar terms from the iteration 
part of the quasipotential. Consequently, to find the necessary contribution 
of order of a^, we must use first approximate equality in (14). Taking the 
difference 



1 



1 



2^2 (po - + ^e) 2m2(po + ie) ^ 4mi(po + ie)' 

let represent the quantity 1/D^(p) in the form: 

1 1 {f + w^) 



(16) 



+ 



(p) 2m2 (po + «e) 4m2 (po + «e) 



(17) 



Second addendum of (17) is of higher order on mi/m2 in comparison with the 
first. But it leads to the necessary order correction on the other parameter 
a. Using the splitting (17), in the sum (15), let extract the terms, which 
generate the correction 0{a^) and 0{a^\na) in the energy spectrum. We 
may write them in the following manner: 



I if + w') 

(p — j/y + w 



2niS{pQ) _ 27ri(5(po - po) _ 2m6{po] 
{po + ieY (po + (Po + if) 

— 27ri5(po) 27ri5(pQ— Po) 2niS{p[ 



+ 



(18) 



(po + ie)2 (po + ie)2 
•2 r 27ri5(p[,) 27ri5(p[, 



+ 



+ 



(Po + «e)^ 
Po) 27ri(5(po 



+ 



(po + ie)2 (po + «e)2 (Po + ie)2_ 

It is evident from three-photon interaction amplitude of the type (6), that 
the parts of (18) give the necessary order corrections on a in the studied fine 
structure intervals. The same order corrections 0(mi/m2), as well as (18), 
will arise from the quasipotential terms containing the functions fii{p,p'), 
when we use the second approximation of (14) for muonic denominators. To 
do definite conclusion about the order of appearing terms in energy spectrum, 
which are determined by these quasipotential addenda, let transform them 
for greater simplification. Let consider for definiteness massless terms in the 
function fii{p,p'), proportional to p^,p'^,pp': 



3p^ 



+ 



DM D^p) D^{p-p') D^-p) 



+ 



(19) 



6 



— 6pp' 



+3/ 
1 



D 



+ 



1 

1 



+ 



1 



1 



+ 



1 



D^,{p'-p) Df,{p-p') 
1 1 



1 



3p 
2m 



D,{-p') D,{p') D^i-p) D^ip) D,{p'-p) D,{p-p' 

- [-27rz5(po) + 27ri5{po - p'o)] + [-27ii5{p'o) + 27ri5{po - p'o)] - 
6pp' 



2mo 



-2TTi5{po 



2m2 

27Tt6{p'o) + 2ni6{po - p'^)] . 



The transformation of other terms in functions fn {p, p') may be carried out 
by analogy. The next period of calculation consists in the integration over 
four-momenta in expressions (18)-(19). The typical two-loop integral, that 
results on this way has the following structure |T6|: 



d'^pd'^p' 



GiiPo,Po,mi)P{p,f. 



w] 



-(27r)8 - t(;2 + ^g) [(p _ p')2 + ^g] (p2 _ ^2 + ie) D ^{p') D ^{p) ' 

(20) 

where Gj(pQ, ''^i) contains one 5 - function, and P{pi,p,w) is a polynom. 
To calculate fundamental integrals (20) we have used Feynman parameteri- 
zation in order to combine the denominators of the particle propagators, and 
the symmetry properties of the integral with the replacement p -v^ p' . There 
are the next set of functions Gi{p'Q,pQ,mi), which appear in this paper: 

2Txi5{pQ — Po)2mi ^ _ 27ii5{po)2mi 



{po + iey 



Go 



ip'o 



G3 = -27riS{po)2mi, 
(21) 

G4 = — 27ri(5(po)2mi, G5 = —2ni6{po — Po)2mi. 
The results of the integrations for Ki (20) have presented in the table |]16|. 



Table of the integrations Ki (20), appearing in the 
muonium fine structure calculations 





p^{pjf) 




ff^iffff) 


w'^ipp) 




21n2- i 


21n2-i 


21n2- i 







ff^ff^ — pff) 


pff {ff^ — pff) 


w'^ijf'^ — ffff) 


2 ->i 
w^p^ 


K2 


1 1 mi 1 

2 2w 32 


1 1 mi 13 

4 2w 32 


5 

32 


2 






ffff 


-f> 




K3 




1 ]„ mi 1 

4 2w 4 


1 ]„ mi 1 

2 2w 8 


1 




1 ]„ mi 1 

2 2in « 


1 ]„ mi 1 
4 2w 4 








In 2 


ln2-i 


In 2 
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Then the contributions, defined by expressions (8-13) and (18) are corre- 
spondingly equah 

AB.^-hzar'^ (22) 

Z 1712 
171^ 1 1 

A52 = (^a)^^(61n2--) (23) 

Wc have introduced in (20) the photon mass w to avoid "infrared" sin- 
gularities. The "infrared" logarithms \nw, containing this photon mass (see 
table of integrals Ki), and appearing at intermediate expressions, are mutu- 
ally cancelled in the corrections ABi, AB2. 

Let consider now the quasipotential addenda, containing the momenta of 
particle relative motion in the initial and final states. We denote them by 
fi and r2 correspondingly. Their consideration leads to modification of fn, 
which acquire the following additional terms: 

A/21 = 10mipV2 + 5pop'r2 + mipr2 + 3popr2, (24) 

A/31 = — mipVi — Spop'ri — lOmipri — Spgpri, (25) 

A/41 = mi(7pVi + 5p'r2 + lOpvi + 5pr2) + Gpop'ri + 5pop'r2 + Sp'^pri + 5popr2, 

(26) 

A/51 = mi (— 5pVi — 10pV2 — 5pri — 7pr2) — 5pop'ri — 8pop'r2 — ^p'^pri — 6pQpr2, 

(27) 

A/ei = mi ( - 1 IpVi - 1 lpV2 - 1 Ipri - 1 lpr2) - 8popVi - 8popV2 -8popr 1 - 8popr2 . 

(28) 

Using again the symmetry properties of appearing integrals under simulta- 
neous variable replacement p <^ p', ri 4^ r2, we obtain cancellation of all 
integrations in (24)- (28). So, the contribution the particle relative motion in 
the fine structure with the accuracy 0(mi/m2) is equal to zero. Thus the 
full value of the calculated correction (ZQ;)^m^/m2 for hydrogen-like system 
S-states is defined as a sum of expressions (22) and (23): 

AB = {Zafi^-^ f61n2 - -V (29) 
^ ^ mim2 V 487 ^ ^ 

Numerical value of obtained contribution (29) for the muonium fine structure 
interval 2^5"! l^Si takes the value 0,271 Hz. We have used in this paper 
the diagrammatic approach to the calculation of the corrections of order 
(ZQ;)^m^/m2. We have made the most complicated part of calculations, 
connected with the two-loop integrations. In order to obtain the total value 
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of necessary order contribution it is important to complete these results by 
similar corrections from one-photon, two-photon interactions as well as by 
the results from the second order perturbation theory [|T^]. It is the aim of 
future studies. 

The work was performed under the financial support of the Russian Foun- 
dation for Fundamental Research (grant no. 96-02-17309). 
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